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Section 1 


1.0 INTRODUCTION 


This program had two objectives: 1) to demonstrate the level of tem- 

poral and fixed pattern noise and the spectral quantum efficiency attainable 
with a Charge Injection Device (CID) imager operating a star tracking mode;' 
and 2) to project the results to a conceptual 400 x 400 CID array for NASA 
star tracker applications. 

Howard of General Electric^) and Stanton et al . of the Jet Propul- 
sion Laboratory(2) have evaluated the potential tracking performance of a 
large CID sensor. The General Electric study used computer simulation and 
measured small CID characteristics to predict accuracy and data rate perfor- 
mance. Both studies concluded that the CID should be capable of meeting the 
application requirements with a minimum of associated equipment. 

The outstanding features of the CID are high quantum efficiency, non- 
destructive readout, and simple interfacing with digital control electronics. 
The X-Y addressability of the CID sensor allows direct access to that portion 
of the array containing the image information. This feature allows operation 
at very low data rates and consequently allows low amplifier noise bandwidth. 
An operator can utilize the NDRO feature to perform repeated readouts of the 
image information or on-chip coding of the pixel quantities, both of which 
can lead to signal enhancement. 

% 

In this program a small (128 x 128) CID array, processed on bulk si 1 i - ' 
con and incorporating transparent upper-level conductors, was evaluated in 
a cooled breadboard camera. All the timing of waveforms and the processing 
of video data was done with a micro-computer. At the same time a system re- 
quirements study was performed to define performance goals and other opera- 
tional requirements. Based upon performance of the large sensor as extrapo- 
lated from the test results, a tentative star tracker system was defined. 
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Section 2 


2.0 PROGRAM HISTORY 

During the course of this program it became necesssary to slightly 
alter the initially proposed technical approach, particularly with regard 
to array fabrication, packaging, and readout methods. 

The proposal to include preamplifier transistors within the array 
package proved burdensome in this experimental test system and was abandoned 
for devices mounted outboard of the array package. While this choice re- 
quired greater scrutiny of electrical interference and increased measured 
output noise (to an unknown degree), it facilitated packaging, handling, 
and testing of devices. More importantly, it permitted direct comparison 
of the imager noise characteristics,, unaffected by preamplifier parameter 
variations . 

The initial readout concept, in which a relatively long time interval 
occurred between the resetting of the output bus and the single signal sample, 
placed a stringent requirement on device junction leakage current, which in 
turn added shot noise to the readout signal. When researchers encountered 
more junction leakage than expected on the earlier processed arrays, they de- 
vised new readout approaches to minimize this effect. Section 3 details a new 
approach in which the array is read out at a higher rate, non-destructi vely, 
and successive digital summations are taken to obtain an average signal. Be- 
cause the processing gain obtained by summation balances the excess thermal 
noise incurred by the increase in bandwidth, system thermal (Johnson) noise 
remains unaffected, while the effect of junction leakage is greatly reduced 
by repeated reset operations at the readout rate. Another benefit of this 
operation mode is that problems associated with low frequency amplification, 
drift, 1/f noise, etc., are equivalently reduced. Subsequently, researchers 
determined and corrected the causes of the high rate of junction leakage, how- 
ever, the summation- and-averaging readout approach was retained for its oper- 
ational ease and other peripheral benefits. 

The third and last change was the substitution of thin («500A ) poly- 
silicon for metal-oxide as an upper level conductor. Although it was intro- 
duced as a "back-up" approach when processing problems with the metal oxide 
were encountered, the polysilicon proved superior to the metal oxide in several 
significant ways. These include ease of fabrication, compatibility with present 
semiconductor processing, and superior spectral performance (an unaccounted for 
result at this writing). The higher resistance of the thin conductors is not 
a performance factor in the subject application because of the low readout 
rates and the fact that readout occurs on the orthagonal (high conductivity) 
lines. To ensure safety and good step coverage, apportion of the conductor 
is maintained at the conventional thickness («5000A): 

In summary, three variations of the initial system concept were intro- 
duced during the course of the program. The first variation, a packaging of 
devices, was done for expediency, and has no relevance to the final system 
concept. The second and third variations are design improvements and Section 6 
includes them in the final system definition. 
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Section 3 

3.0 SYSTEM DESCRIPTION: CID ARRAY DESIGN 

3.0.1 Cell Layout 

Figure 3-1 shows the 30 x 30 pm cell layout of the 128 x 128 sensor. 

The thick-thin oxide boundary is dashed, the lower conductor (N-doped-poly- 
silicon) crosshatched, and the upper level (N-doped polysilicon) conductor 
is shown in solid outline. The lower level poly-thin oxide capacitor has 
an area of 1.36 x 10 -o cnr fi and the upper level poly-thin oxide capacitor 
has an area of 1.81 x 10 cm . Since the upper level capacitor has a thicker 
dielectric than the lower level capacitor, the charge storage capacity of 
each electrode is approximately the same, 1.9 x 10 carriers. Depletion 
capacitance loading limits the maximum output charge to 1.5 x 10 carriers. 



Figure 3-1. Cell Layout 


3.0.2 Array Layout 

The array in Figure 3-2 contains 128 rows and 128 columns, selected 
by groups of four on both the horizontal and vertical axes. Since the de- 
sign of this array was for sequential readout of the image subblocks, selec- 
tion is by means of scanning registers. In operation, a logical "one" would 
be entered into each scanning register and shifted, as required, to select 
the desired subgroup. The column drive lines, through E^, would then 
be driven to obtain four parallel outputs. ' f 
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In addition to the 128 active rows, another row is available to pro- 
vide differential cancellation of column drive interference.. This compen- 
sation row is selected for every row address and is cleared when every row 
is cleared. 

Overall chip size is approximately 5x5 mm. Figure 3-3 shows a photo- 
micrograph of the array. 

3.0.3 Array Parameters 

Table 3-1 shows the calculated array electrical parameters for the 
128 x 128 format. 
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Figure 3-3. Array Photomicrograph 
Table 3-1 

ELECTRICAL PARAMETERS 


Cell Storage Capacity 

1.9 x 10 carriers 

Saturation Output Charge 

1.5 x 10 carriers 
(.24 pC) 

Distributed Row Capacitance (C RQW ) 

6 pF 

Sense Line Capacitance (C<0 

12 pF 

Distributed Row Resistance 

3.9 K ft 

Row Select Transistor Resistance 

4 Kft 


3.1 READOUT METHOD 

The readout technique for the track mode noise evaluation is the "double 
read" method described by Howard (1). As Figure 3-4 illustrates, a 4 x 4 
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| “background" signal plus 

“background" 

inject inject inject 

TIME 

Figure 3-4. Track Mode Background Rejection Technique 

pixel group is read out twice, non-destructively, and the first (background 
group subtracted from the second group to form the net signal. This 
scheme effectively cancels the unwanted background signal (pattern noise) 
which repeats for each reading; however, a penalty is incurred because sig- 
nal integration time is halved and Johnson noise increases. These effects 
will be discussed in more detail below. There are other approaches for 
background removal that could be implemented with less noise penalty. These, 
however, involve system considerations that could not be adaquately evaluated 
under this program. It was decided to conduct the evaluation in this envi- 
ronment which unequivocally addresses the pattern noise concern even though 
the resulting signal-to-thermal noise ratio is degraded. 

The noise performance to be reported here was obtained using the con- 
ventional sequential (block code) readout of the individual pixels. Hadamard 
and other transform readout methods were not employed. Transform readout 
remains an option for future designs where system operating characteristics 
are better defined and/or the need for a further reduction in system thermal 
noise is obvious. 

3.2 CIRCUIT DESCRIPTION 


Figure 3-5 shows a block diagram of the complete control and signal 
processing system. All control and timing waveforms were computer gener- 
ated; and the same computer, using a separate program, processed the stored 
video data, calculated noise levels, etc. 
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Figure 3-5. System Block Diagram 
3.2.1 Preamplifier Circuit 

Since amplifier noise is a significant factor in the overall noise 
inventory, the design of this function is critical to system performance. 
Particularly for the star tracker application, low frequency (1/f) noise 
is a concern. Therefore, on-chip MOS devices were discarded for off-chip 
J-FET's as input devices. The latter exhibit significantly lower noise 
at low frequencies, although their use increases input capacitance. The 
particular choice (2n6451) depended on the availability of these devices 
in chip form for mounting in the array package, and on the drain-to-gate 
leakage. Later changes in packaging and operating strategies subsequently 
disposed of both requirements; however, the program schedule did not permit 
circuit redesign. 

A schematic diagram of the preamplifier circuit is shown in Figure 3-6. 
The input transistors operate as common-drain (source-follower) amplifiers. 

The outputs of these devices are dc coupled to bipolar differential -pairs 
which form the difference taking function between the compensation line 
and each of the four signal lines. A final stage operational amplifier 
provided double-ended to single-ended conversion of the amplified difference 
signals. Direct coupling was used throughout to avoid the problems associated 
with the use of coupling capacitors in the low impedance, low frequency 
environment. 

In addition to the design of the basic preamplifier circuitry, great 
care was exercised in the design of peripheral circuitry, such as power 
supplies and biasing networks, to avoid contamination of the signal with 
excess thermal and induced noise from these sources. Finally, care was 
taken in the use of electronic instrumentation, including the data acquisi- 
tion and computational circuitry, to minimize the noise contribution of 
the 60 Hz power lines. 
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Figure 3-6. Preamplifier Schematic 

3.2.2 Post Amplifier 

The signal level after preamplification is high enough so that conven- 
tional data acquisition circuitry can function without significant contri- 
bution to the overall system noise. As Figure 3-7 shows, this circuitry 
consists of variable gain stage, used primarily to balance the gain of the 
individual channels, and two sample-hold circuits. In effect, the first 
S/H performs the restore function, while the second device samples the video 
signal and- holds it for multiplexing into the A/D converter (not shown). 

3.3 TIMING WAVEFORM GENERATION 


The timing waveforms needed to drive the array in the "double read" 
method were sufficiently slow so that an assembly language program written 
for a 4 MHz Z-80 microcomputer could deliver all necessary signals at the 
proper times. A block diagram of the microcomputer. Figure 3-8, shows that 
the 11 digital signals required for array operation are sent to the driver 
board by outputting an 8-bit number to 1 of 3 latched parallel output ports 
of the Z-80 system. Thus, 11 of the 24 available bits are used for array 
timing. The drivers and voltage regulators on the driver board converted 
these TTL level signals to the appropriate MOS levels. Some of the other 
available bits are used directly for functions such as restore, sample and 
hold, A/D conversion, and sync pulses for the oscilloscope. 

Figure 3-9 shows the general timing waveforms used in the "double 
read" method. The four columns in the 4x4 block are read sequentially, 
while the four rows are processed in parallel. After injection, each column 
is read for background and then 50 ms later it is read again for background 
plus signal. The first reading is subtracted from the second so that each 
column sees only 50 ms of signal integration time. 
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Figure 3-7. Post Amplifier/DC Restore/Sample 

■ Figure 3-10 shows the injection timing in more detail. After reading 
of a block is complete, the vertical phase lines 4>1V, 4>2V and the horizontal 
phase lines <j>lH, <()2H are run out to clear the scanners. Then the array 
is injected by driving RVD and CVD simultaneously. After a recovery period, 
the vertical and horizontal inputs (Vin, Hin) are applied to the vertical 
and horizontal scanners and the phase lines are run until the proper block 
is selected. The array is then ready for the read sequence. 

Figure 3-11 shows the read sequence in detail. In actuality, every 
time the computer reads a column, as per Figure 3-11, it reads that column 
64 times and averages the readings. This method minimizes the leakage cur- 
rent component of noise in comparison to the other noise sources. It 
also allows the amplifier to be used beyond its 1/f noise region. This 
will be explained further in Section 4. 

For each of the 64 readings of a column, pulsing IG negative resets 
the output lines; and any dc offset from the postamplifier is taken out by 
the restore pulse. Next, the appropriate E-line is driven, and after 4 ampli- 
fier time constants ( 80ms) the output level is sampled and held simulta. - 
neously on all four rows. While the reading is being taken, and computer 
and A/D converter process the previous reading. In addition to the 64 readings. 
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Figure 3-8. Microcomputer Block Diagram 
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Figure 3-9. Overall Timing Diagram 
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STAR TRACKER TIMING DIAGRAM 
INJECTION SEQUENCE 
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Figure 3-10. 

Timing Diagram-Injection Sequence 
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Figure 3-11. Timing Diagram-Read Sequence 

one initial "dummy" reading is used but not processed in order to avoid 
start-up transients. After one column is read 64 times, the average reading 
for each of the four rows is stored in computer memory and the next column 
is read. 


A/D READ 

AND SUM ‘ ^ 
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3.4 DATA COLLECTION AND PROCESSING 


The Z-80 microcomputer collects and processes data and provides the 
timing waveforms for the star tracker. A 12-bit A/D converter module with 
an 8-channel differential -input multiplexer, buffer amplifier, and sample- 
and-hold circuit is used at its maximum throughput rate of 100 kHz. After 
the' sample is taken and held for the four rows of each reading, the computer 
reads the four rows by sequentially updating the multiplexer channel and 
making conversions. 

The computer stores an 18-bit sum of the sixty-four 12-bit readings 
for each pixel in the 4x4 group, for both the background and background- 
plus-signal. The 18-bit difference is then stored in the proper location 
in memory. We have determined that the converter quantization noise is in- 
significant as long as the RMS value of the noise to be measured is greater 
than, or equal to, the value of one least significant bit of the converter. 

For our setup, this LSB corresponded to 200 carriers. Another advantage of 
averaging 64 reading is that the noise on the average is (64)1/ 2 or 8 times 
less than the noise on the individual readings. As a result, we can obtain 
the precision of a 15-bit converter while using only 12 bits. The difference 
taking adds a factor of -J2 to the noise, since the noise adds incoherently. 
Thus one can measure noise levels as low as 35 carriers RMS with no signi- 
ficant quantization noise contribution. 

Several versions of the computer program were written t.o perform the 
various tests. One version read out the whole array, a 4 x 4 block at a 
time, and stored the result in memory so that it could be displayed on a 
television monitor using the direct memory access display board shown in 
Figure 3-8. It was possible to simultaneously store and display two pictures, 
each containing an arbitrary 8-bit section of the 18-bit difference computed 
for each pixel. For example, one could show the most significant eight bits 
when a star was imaged or the next eight bits when the lens were closed, in 
order to see the background fixed pattern noise. This allowed a visualization 
of the tremendous dynamic range provided by the star tracker. 

Another version of the program was written to measure the temporal 
noise level. This program chose the center block of the imager and read 
the block 100 times, storing the result each time for each of the 16 pixels 
in the block. A separate program, written in BASIC, calculated the mean 
and standard deviation (or RMS noise) for each pixel. Since the sample 
variance decreases linearly with the number of readings, the RMS error on 
these noise readings is 5%. Finally the program takes the average mean 
and standard deviation over the 16 pixels, thus reducing the error. These 
are the numbers that the noise analysis section quotes. 

For the total. noise measurement (fixed pattern and temporal noise) 
a slight variation was made in the temporal noise program. In this case 
100 different blocks around the array were read and" stored, and the same 
BASIC program was used to analyze the data. 
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Figure 3-12 shows a typical printout of total and temporal noise, 
respectively. 
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Figure 3-12. Typical Noise Data Printout 
3.5 TEST FIXTURE DESCRIPTION 


All tests were carried out and data taken with the imager and preampli- 
fier mounted in a demountable liquid nitrogen dewar (Infrared Laboratories 
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Inc., Model HD-2). Temperature of the dewar work surface was maintained by 
a controlled flow of chilled dry nitrogen into the dewar in place of the normal 
procedure of using liquid nitrogen directly. The dry nitrogen was cooled by 
heat exchange from a supply of liquid nitrogen in an open dewar, and the array 
temperature was monitored by a digital thermometer via a thermocouple attached 
to the imager package. Drive circuitry and post amplifier circuits were mounted 
on the dewar assembly and wired to the internal circuitry through hermetic con- 
nectors. An optical window in the dewar wall permitted the imaging of a light 
spot for verification of proper system operation. Figure 3-13 is a photograph 
of the complete system. 
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Section 4 


4.0 NOISE PERFORMANCE 


The circuitry used t.o select and provide readout of CID image sensors 
contains a number of Johnson noise, sources. The distributed resistance of 
the array lines used for signal, sensing, the line selection switch, and the 
first preamplifier stage each contribute temporal noise to the video signal. 

In addition, capacitor reset noise (KTC noise), shot noise in the dark cur- 
rent, and/or junction leakage current in the MOS line select multiplexers 
can be significant under certain conditions. 

If a switch is used to set the voltage across a capacitor, thermal 
noise in the resistance of the switch results in an uncertainty in the final 
capacitor voltage. The magnitude of this uncertainty [3] is 

V n « (KT/C) 1/2 (5) 


or 


Q n = (KTC) 1/2 (6) 

where K=Boltzmann 1 s constant = 1.38 x 10~ 23 W-s/°K and T = temperature in 
degrees Kelvin. It is possible to reference the charge signal to the input 
capacitor voltage after reset has been completed. This technique, called 
correlated double sampling results in the substitution of KTC noise 
on a clamping capacitor for KTC noise on the array output capacitance. The 
level of KTC noise referred to the array can be made arbitrarily small, 
however, since gain can be used between the array output and the clamp 
capacitor. Voltage noise at the input of the preamplifier results in an 
equivalent input charge that is directly proportional to the array output 
capacitance (q = cv). Theoretical preamplifier noise levels of a few hundred 
carriers result from array output capacitance levels in the 10-pF region. 

Under low video rate readout conditions, Johnson noise can be minimized 
by restricting the noise bandwidth of the video amplifier. Shot noise orig- . 
inating in array dark current and line select multiplexer junction leakage 
can be limiting under these conditions. Reductions in array operating tem- 
perature can control these thermally generated currents and, consequently, 
the resultant shot noise. 

4.1 THEORETICAL TEMPORAL NOISE LEVELS 


The theoretical levels of the significant noise sources in CID image 
sensors have been investigated and experimentally verified in a detailed 
noise study. (4) The significant temporal noise sources encountered in 
an array operating in a star-tracking mode are defined below: 

a) Distributed Row Resistance 
N DR = (4 KTBR r / 3) 1/2 C R /q 
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where 


N DR 

K 

T 

B 

R R 

C R 

q 


rms noise equivalent carriers 

Boltzmann's Constant 
absolute temperature 
noise bandwidth 
row resistance 
row capacitance 
electronic charge 


where 


b) Row Select Switch 

N rs = (4KTBR s ) 1/2 C R /q 

N rs = rms noise equivalent carriers 

R c = Row Select Switch resistance 

s » 


The channel resistance of an MOS. switch, R , is given by: 


where 


R c = l/(pC 0 (W/L)(Vgs-V t )) 


P 

C o 

W 

L 



= carrier mobility 
= oxide capacitance 
= channel width 
= channel length 
= gate-source voltage 
= threshold voltage 


Parasitic source and drain resistance resulting from the non-zero sheet re- 
sistance of the transistor source and drain region must be added to the c 
channel resistance, Rc, to obtain total selection switch resistance, Rs. 


c) KTC Noise 

N KTC = ( ktc f) 1/2 ^ = 401 c p 1/2 


where 



= rms noise carriers 
= capacitance. Farads 
= capacitance, pico Farads 


The degree to which KTC noise is surpressed by correlated double sampling 
is determined by the degree to which the first sample is a measure of the 
KTC offset voltage. For a video system whose bandwidth is limited by a s 
single time constant, the magnitude of the sample voltage is: 


V , = V . , (l-e" T/RC ) 

sample signal 
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If the sample interval. 


^ample^signal 
factor of 50. 


0.98. 


T, is made to equal 4 time constants, then 
In this case, KTC noise would be surpressed by a 


d) Amplifier Noise 

The equivalent noise resistance of a junction PET is: 
R a = 2/(3 ?m> 

where 


q = transistor transconductance 
The amplifier noise is given by: 

N a = (4KTBR a ) 1/2 C t /q 
where N = rms noise carriers 

ci 

R = Transistor noise resistance 
a 

C t = Total input capacitance 

e) Junction Leakage Current Shot Noise 

Nj - dj^M ) 172 

where N.-, = RMS noise carriers 

J i 

I., = junction leakage current 
J i 

t = sample time interval 

f) Dark Current Shot Noise 

N D = ^ 

where Np = RMS noise carriers 

Ip = dark current per pixel 
t^ = integration time interval 

Table 4-1 lists the theoretical noise levels for the FPP-128 array 
operated in this star tracker evaluation camera. 
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Table 4-1 

THEORETICAL NOISE LEVEL 


Source 

Spot Noise 

Effective 
Capacitance - 

Noise Level 


(nV/Hz 1/2 ) 

(pf) 

1/2 

(carriers /Hz ' 

Row Resistance (5Kfi) 

5.2 

6.0 

0.20 

Row Select Switch (3K&) 

7.0 

6.0 

0.26 

Disconnect Switch (2.4K&) 

6.3 

10.4 

0.41 

J-FET Preamplifier 

1.1 

31.0 

0.21 * 

Differential Amplifier 

1.9 

31.0 

0.37* 

(per input) 






Root Mean 

Square 0.68 


•^Theoretical RMS Amplifier Noise (0.43) compares with measured value of 0.57) 


The video amplifiers operated with a bandwidth limiting time constant- of 
20 microseconds. This gives a noise bandwidth of B = ( tt/2 ) ( 1 /2ttRC ) = 12.5 kHz, 
Correlated double sampling results in the addition of amplifier noise power 
for each of the two' samples, effectively doubling the noise bandwidth to 
25 KHz. Each of the two correlated sample intervals were set at 4 time • 
constants (80 ps) to suppress KTC noise by a factor of 50. Leakage current 
shot noise was neglible over this sample time interval. 

4.1.1 Low Frequency Operation 

The dominant temporal noise at high video rates is from Johnson noise 
sources (resistors, switches, field effect transistors). As the operating 
frequency is reduced a reduction in Johnson noise, proportional to the square 
root of bandwidth, can be realized. At lower rates, however, low frequency 
excess noise and leakage current begin to interfere with the sensing opera- 
tion. It would appear that these low frequency npise effects would dominate 
below some video rate and limit the improvement which this approach could 
achieve. This is not the case, however. There is a strategy that allows 
these low frequency effects to be circumvented and the effective noise band- 
width to be reduced with, at this writing, no known limit. 

The strategy is to choose the imager operating frequency for lowest 
noise, somewhat above the point where low frequency effects become signif- 
icant, and then use multiple non-destructive readout operations and external 
signal summation to further reduce the effective noise bandwidth. Because 
signals sum coherently and uncorrelated noise incoherently, the effective 
noise is reduced in proportion to the square root of the number of non- 
destructive readout operations. This strategy is effective in reducing 
all temporally uncorrelated noise, including Johnson noise, KTC noise, charge 
transfer noise, and leakage current shot noise. 
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4.2 FIXED PATTERN NOISE 


Solid-state image sensors can exhibit a fixed non-uniform spatial 
background in the reproduced image. The major sources of fixed pattern 
noise in CID image sensors are transistor switching interference, array 
photolithographic variations, bias charge variations, and silicon material 
defects . 

One particular silicon material defect that could affect tracking 
accuracy is resistivity striations in the single crystals grown from the 
molten, state, (5 ) These striations could lead to lateral movement of 
photon-generated charge between the location from which it generated and 
the pixel storage region. The use of more uniformly doped epitaxially grown 
silicon or transmutation-doped silicon should avoid this potential problem, 
which is common to all silicon image sensing devices.- Non-uniform coupling 
of the MOS transistor scanner output voltage to the video signal results 
in a component of fixed pattern noise that repeats from scan to scan. Varia- 
tions in row-to-column crossover capacitance arising'from either insulator 
thickness or photolithographic variations cause a two dimensional component 
of fixed pattern noise. Variations in the bias charge from site-to-site, 
caused by differences in storage capacitance or threshold voltage, also 
result in a two-dimensional component of fixed pattern noise. Differential 
sensing and signal processing can minimize these effects. 

Dark current non-uniformity can be an important source of pattern noise, 
particularly at room temperature. The inherently low dark current performance 
of CID imagers is an advantage under these conditions. 

Site-to-site variations in sensitivity will result in the appearance 
of FPN which is more pronounced in the highlight portion of the image. 

This effect could be the result of patterning variations in any opaque or 
semi -opaque layers on the array surface. In some instances, this might 
result from gain differences of amplifier channels that are subsequently 
multiplexed onto a single video line. The latter cause can be corrected 
by suitable adjustment, and is not considered a fundamental FPN source. 
Sensitivity variations of the FPP-128 array were measured by first finding 
the mean signal for 25 readings at 3/4 full-scale illumination for each 
of 50 sites along a row. Similar values, measured for no illumination, 
were subtracted from these values. The 49 differences of these signal-minus- 
background values were then averaged to determine the pixel-to-pixel sen- 
sitivity variation. Measured sensitivity variations of 1/2 percent are 
typical for these imagers. Thus, the FPN contribution of sensitivity vari- 
ations is relatively small. 

4.3 MEASURED 'RESULTS 


Temporal noise was measured as a function of temperature with the 
imager operating in the track mode. In this operational mode a 4 x 4 sub- 
group of pixels is read out twice each frame. The difference between these 
two readouts is taken to cancel fixed pattern noise. An integration time 
of 50 ms results from this procedure, with the imager operating at 10 frames 
per second. A single 4x4 sub-group was read out 100 times at each temperature 
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setting and the standard deviation of the readings for each pixel was cal- 
culated. The average of the 16 resultant standard deviations is plotted 
as the temporal noise level at each temperature. 

Total noise was measured by reading a single frame for each of 100 
different 4x4 sub-groups distributed over the image sensor. The standard 
deviation was calculated for each pixel of the 4x4 sub-group. The average 
of the 16 resultant standard deviations is plotted as the total noise level 
at each temperature. 

Figure 4-1 shows the temporal and total noise measured on three arrays 
from wafer #9-7. This wafer exhibited high dark current which is respon- 
sible for the high total noise observed. The spatial variation in dark 
current gives rise to a pattern noise which varies with temperature in the 
same manneg as the dark current itself, which decreases by a factor of two 
for each 8 U C reduction in temperature. Shot noise associated with this 
dark current is proportional to the square root of the quantity of dark 
charge collected and consequently is reduced by a factor of 2 for every 
16 C reduction in temperature. This effect can -be seen in the temporal 
noise data for imager #9-7-15. Figure 4-2 represents results obtained with 
arrays from a wafer with the proper dark current level. Dark current shot 
noise ig lower than the Johnson noise level of 40 carriers for temperatures 
below 5 U C. The spatial variation of dark current, the dominant fixed pattern 
noise at room temperature, is lower than other fixed pattern noise components 
at 0°C. The increase in total noise at very low temperature is a fixed 
pattern effect caused by the scanning registers. It is believed that certain 
scanner circuit nodes, which do not have their voltage levels well-controlled, 
caused the effect. This effect can be avoided with a different scanner 
design or with a digital decoder for row and column selection. 
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Figure 4-1. Measured Noise-Wafer #7 
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Figure 4-2. Measured Noise-Wafer #12 
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Section 5 


5.0 SPECTRAL RESPONSE/QUANTUM YIELD 

In the following presentation, the term "quantum efficiency" refers 
to the ratio of charge collected to incident optical photons (see Figure 5-1). 
In sensing this stored charge in a CID imager, a loss is incurred due to 
the charge lost in restoring the depletion region under the storage electrode 
after injection. For the subject device, the factor was 0.82. Here the 
term "quantum yield" will describe the ratio of detected charge to incident 
photons. As indicated in Figure 5-1, any additional gain losses are con- 
sidered part of the overall system gain. 

CHARGE CHARGE SIGNAL 



QUANTUM OVERALL SYSTEM 

YIELD GAIN 


Figure 5-1. Signal Flow Diagram 

The data that follow are presented in. terms of quantum yield in order 
that it be referenced to the same point in the signal chain as the noise 
data presented above. 

5.1 MEASURED RESULTS 

The quantum yield of six imagers having thin polysilicon upper electrodes 
was measured over the spectral range of 0.41 to 1.05 pm at room temperature. 
Each imager .was also cooled while measuring quantum yield at .8 pm to evaluate 
any temperature effects on imager spectral performance. 

Measurements were made by placing the imager in an evacuated dewar and 
illuminating the array through a dewar window with light from a constant energy 
glass prism momochromator (Perkin Elmer Model 98). At the imager location a 
silicon photo. diode (GE #E07) measured incident light power per unit area at 
each functioning wavelength. An NBS calibration (Test #212736) determined its 
absolute spectral response. 

During charge sensing, electronics located outs-ide the dewar maintained 
all row and column electrode voltages at -10V. This permitted simultaneous 
collection of photo-generated charge at each sensing site in the array. The 
value of the resulting- substrate current was used to calculate quantum yield, 
which is defined as the ratio of the number of electrons read from the imager 
to the number of incident photons. 
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Results of the quantum yield measurements for imagers from wafer #9-7 
and #9-12, appear in Figures 5-2 and 5-3 respectively. The average quantum 
yield of the six devices is 24% at .41 urn and above 45% from .6 to .95 pm. 



WAVELENGTH (MICRONS) 


Figure 5-2. Quantum Yield vs Wavelength: Wafer #7 

Quantum yield measurements were also made for imagers having ATO upper 
electrodes. At corresponding wavelengths, the ATO imagers have a 5-25% lower 
QY than the thin polysilicon imagers. 

Failure of the transparent electrode material to produce higher quantum 
efficiencies cannot be accounted for, but an excessive charge recombination 
rate at the silicon surface could be responsible. 

5.2 TEMPERATURE EFFECTS 

Cooling the imagers to -20°C produced no measurable change in quantum 
yield at .8 pm. However, studies of the silicon energy absorption curve 
at 77°K (-196°) by Dash and Newman L 6 ) show a temperature coefficient of the 
bandgap of about -4 x 10“ 4 eV/°C. These data suggest that the quantum yield 
curves in Figures 5-2 and 5-3 would show a significant leftward shift if 
measured at 77°K (-196°C). The shift would be more pronounced at longer 
wavelengths, decreasing QY by 30% at 1.0 pm; but the shift would have little 
effect at wavelengths shorter than .85 pm. Over the temperature range of 
'-20 to +25°C; however, this effect is negligible. 
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WAVELENGTH (MICRONS) 

Figure 5-3. Quantum Yield vs Wavelength: Wafer #12 
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6.0 PROJECTED PERFORMANCE OF LARGE ARRAY 


Projecting the performance of .the 128 x 128 test array to an array of 
order of 400 x 400 involves more than simply proportionally scaling the 
test array parameters to those of a larger device, because the. larger device 
would differ significantly from the test array fn both design and. level 
of technology.. In addition, changes in readout methods and circuitry would 
affect final performance. " . 

While at this time the actual configuration of the 400 x 400 array's 
selection techniques, tracking sub-group size, and other system-related 
details is uncertain, some changes in design are clearly indicated. The 
unneeded "disconnect" transistor, a major noise source in the test device, 
would be eliminated and capacitive coupling between the signal lines and 
compensation line would be reduced. A symmetrical cell layout is also de- 
sirable to minimize signal irregularities caused by photomask mis-registra- 
tions. Finally, improvements in processing technology would reduce the 
sensing site size and, hence, overall chip size. A 20pm square sensing 
site size, for instance, would result in a 8 x 8 mm sensing area with the 
overall chip size being only slightly larger. This size chip could be pro- 
cessed with high confidence and with excellent projected performance as 
will be detailed below. It should be noted that this choice of array size, 
in the absence of specific system requirements, is somewhat arbitrary. 

From a system perspective, the preference is for a large number of sensing 
sites. In a multi -array system, there will be a trade-off decision among 
the number of arrays, array cost, and complexity of manufacture. The results 
derived here can be considered a reference point for any final system decisions. 

6.1 NOISE 


As discussed in Section 4, each on-chip noise- source must be considered 
in relation to the effective capacitance upon which it operates in order 
to assess its contribution to total noise charge. It is necessary, therefore, 
to calculate both the resistance and capacitance of the new design in order 
to establish new theoretical noise values. Using this approach, the noise 
contributions of the four major noise sources were computed and tabulated 
in Table 6-1. It was assumed that double sampling and cooling significantly 
reduced KTC noise and leakage current shot noise. Column 1- of Table 6.1 
lists the baseline values for the 128 x 128 test array. The noise contributions 
of the on-chip sources are nominal values extrapolated from measurements 
of appropriate test devices. The' test fixture directly measured preamplifier 
noise and a total of approximately 10 noise equivalent carriers (NEC) per 
channel of the input differential preamplifier was computed. This count 
yields a total measured value of about 28 NEC, as measured in the test 
fixture, because a noise multiplier of 2(2) ; is incurred in. the various 
differencing operations inherent in the readout method used. (Actual measured 
values, as reported in Section 4, were on the order of 36 NEC.) 
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Table 6.1 

THEORETICAL NOISE/CHANNEL 

FOR 400 x 400 ARRAY FOR DIFFERENT DESIGN CONDITIONS 


Noise Source 

128 x 128 . 
Basel ine 
(NEC//Hz) 

400 x 400 
(NECVHz) 

400 x 400 
(rows strapped) 
(NEC/^/Rz) 

400 x 400 
(N-Channel 
process) 
(NEC//Hz) 

Row 

Resistance 

0.20 

0.75 

0.0 

0.0 

Row Select 
Transistor 

0.25 ' 

0.30 

0.4 

0.22 

Disconnect 

Transistor 

0.44 

0.0 

0.0 

0.0 

Pre-Amp! if ier 

0.57 

0.21 

0.21 

0.21 

Total Spot 
Noise 

0.79 

0.84 

0.45 

0.30 

Total Noise 
(0 B.W. = 150 Hz) 

9.7 NEC 

10.2 NEC 

4 

5*5 NEC 

3.7 NEC 


Column 2 of Table 6.1 gives the equivalent values for the 400 x 400 
array. These numbers reflect not only the effects of scaling, but also 
the effects of actions which can correct deficiencies in the test structure. 
Detailed examination of Column '2 reveals that the noise contribution of 
the distributed row resistance has increased significantly due to an increase 
in its effective resistance and capacitance.. The larger row capacitance 
also raised the noise charge of the row select transistor. The so-called 
"disconnect transistor,” an unneeded feature, is presumed eliminated, and 
the preamplifier noise reduced, the latter by the improved circuit design 
and the signal-averaging readout method. This method permits a higher 
sampling rate. The net result of about 11 NEC per channel is only slightly 
different from that of the test device, an encouraging result. 

The dominance of the row resistance component of this noise,, however, 
leads one to consider the possibility of reducing that effect. The application 
of a conductive "strap" to the offending conductor can virtually eliminate 
this resistance. Column 3 of Table 6.1 gives the result of this action. 

Because of layout constraints, row capacitance must be increased slightly, 
and this change adversely affects the noise charge associated with the row 
select transistor. In addition, some obscuration of the incident radiation 
would occur (about 12.5%). In view of the trade-off- considerations involved, 
this feature might best be considered a design option. , - 


; Ji- i O’ . Xj'iAjt} Wv-; 
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Finally, Column 4 reflects the further noise reduction possible with 
the addition of N-channel processing, currently under development at Gen- 
eral Electric. The viability of this development depends on program timing 
for any projected application. Because of program timing and the realiza- 
tion that an attempt to physically attain the indicated noise level would 
require a re-appraisal of many noise sources now considered negligible. 
Column 4 results are discounted in the system evaluation to follow. 
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Section 7 


7.0 SYSTEM STUDY RESULTS 


The application studies covered in this section of the report were 
supplementary efforts to cover the CID's performance when used in stellar 
sensor applications. The requirements for several stellar sensor applications 
were considered. However, the Space Telescope FGS was selected as a basis 
for identifying the range of sensor appplication since it provides the best 
example of the three primary modes of. operation, namely, acquisition, tracking 
and astrometry. General Electric has already -done significant work on the 
star tracker appl ication. (1 ) This report also considers the use of the 
CID sensor, with an emphasis on the acquistion mode and further evaluation 
of the track and astrometry modes. - 


The net results of the application studies shows the CID continues 
to hold excellent promise for star tracker applications because of its unique 
characteristics. The CID baseline design characteristics are shown in 
Table 7-1.. 


Table 7-1 

CID BASE LINE DESIGN 


1. Chip Size' 

400 x 400 

2. Pixel Subtense (arc-s) 

0.5 

3. Pixel Size (pm) 

20 x 20 

‘4. Total Active Area (mm) 

8.0 x 8.0 

5. Site Storage Capacity (carriers) 

10 6 

6. Effective Capacity (carriers) 

8 x 10 5 

7. Estimated NEC (carriers) 

25 


The results of the requirements analysis task show that the specific 
readout technique used, and its impact on the detail chip design, should 
be selected through a specific set of requirements in order to demonstrate 
the performance advantages of the CID in general and a selected design in 
particular. The different readout techniques can be incorporated into the 
chip design with high confidence once the device is adequately characterized. 

7.1 REQUIREMENTS SUMMARY 

Representative stellar sensor requirements were examined for five 
potential applications for the CID sensor. These applications included 
the Fine Guidance Sensor (FGS) for the Space Telescope, the Fixed Head Stellar 
Tracker (FHST), and a representative booster inertial navigation system 
update case. The results of this review are shown in Table 7-2 where the 
parameters of interest are tabulated for comparison purposes. It is cautioned 
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Table 7-2 

SUMMARY COMPARISON OF REPRESENTATIVE STAR TRACKER REQUIREMENTS 

NOISE MODELS 



, 

FGS/S.T. 

FHST 

STELLAR 

HAADS 

Booster 


Key Parameters 

Aca. 

Track 

Astr. 

Acq. 

Track 

Acq. /Track 


INS Update 

1. 

FOV 

2 

109 min 
(multi chip) 


20 min 2 

Total 8.5° 

O 

CO 

X 

o 

CO 

2° 

8° 

l 

1° x 1° 

r 2 - 

Prob. of Det. % 

85 

- 

- 

95 

- 

- 

95 

- 

t;'3. 

Visual Mag. (M y ) 

7 - 14 

7 - 14 

7 - 21 

2 - 5.7 
Gov 

2 - 5.7 
Gov 



1 - 3 


Eff. Aper. Dia (cm) 

200 



5.8 


- 

2.5 

- 

V'5. 

LOS Rate arc-s 

<2 

<0.2 

<.2 

< 1000 ■ 

< 1000 

- 

6.8 to 
285 

neg. 

,, 6. 

Min. Data Rate 
(Hz) 

- 

20 

1 Meg Bits/s 

- 

- 

3 to 10 

> 1/hr. 

2.6 

f 7. 

Integration (s) 
Times 

>.22 

.025 

,050‘ to 600 

- ^ 

- 

“ 


.04 

8. 

Acq. Time (s) 

60 

- 


10 

- 

0.2 

- 

- 

9. 

£ 

Pos. Acc. (arc-s) 

0.5 

.Ol(lo) 

.01 (la) 

- 

180 (uncal) 
10 (lo) 

1.7 

1.0 

(Bias) 

-1 

* 10. 

Stability (arc-s) 

- 

.0028 

(la) 

.0028 

‘(Iff) 

- 

16 for M =5.7 
7.3 foV M v 

0.1 to 0.7 
M v = 6 

1.7 RMS 


11. 

Mag. Acc. (&) 

- 

- 

1 


- 

- 

12.5 

- 
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that the numbers shown were taken from a variety of sources and have not 
been fully verified in terms of their interpretation or their accuracy. - 
For this reason a number should be considered as identifying representative 
ranges of interest rather than specifications for the application shown. 

For the cases shown, the tightest accuracy and stability specifications 
are for the Fine Guidance Sensor. In practice, the difficulty of meeting 
accuracy requirements is a function of the pixel or chip field of view. 

The numbers can be selected arbitrarily based on either the number of pixels 
per chip or the number of chips used to cover the desired field of view. 

In considering the signal level for a variety of applications, a com- 
bination of several factors must be considered. The principal ones involve 
the effective aperture size, the magnitude of the star to be detected, and 
the integration time available, which is limited by the line-of-sight rate 
and the selected subtense angle of the pixel. A comparison has been made of 
the worst case signal level which might be expected for three of the appli- 
cations shown in Table 7-2. The results are shown in Table 7-3 and indicate 
two limiting cases of interest. The first case shown is that in which the 
•integration time is limited by the pixel subtense angle divided by the line- 
of-sight rate. The most stringent requirement shown would be the High Alti- 
tude Attitude Demonstration System (HAADS) where the signal level is one the 
order of one-tenth of the signal level in the Space Telescope acquisition mode 
because of the field of view and high .spin rate (LOS rate). However, if the 
data frame rate is selected as- the limiting time interval, the signal for the 
FGS acquisition mode is about the same as HAADS. The signal levels shown 
for the astrometry mode can be selected based on the desired integration time 
and hence do not represent a limiting case. These figures indicate that for 
sensitivity as well as accuracy, the .FGS is the most stringent requirement. 
Hence, in the analyses which follow, the FGS for the Space Telescope require- 
ments has been used as a basis for the application potential of the CID. A 
summary of the assumed requirements for each of the^three modes is presented 
in Figure 7-1. 


Table 7-3 

REPRESENTATIVE WORST CASE REQUIREMENTS 


App./Mode 

Eff. 

Ap.Dia 

(cm) 

Pixel 

Sub 

Tense 

(arc-s) 

Visual 

Mag. 

(™ v ) 

Los Rate 
Limited 
t.* 
(s) 

Relative 

Signal 

Strength 

Frame 

Rate 

Limited 

(s) 

Relative 

Signal 

Strength 

Space Tel/Acq. 

200 

.5 

14 

.25 

1 

.025 

0.1 

Space Tel/ 
Track 

200 

.5 

14 

2.5 

10 

.025 

1.0 

Fhst/Track 

5.8 

60 

6 

.07 

.3 

- 

- 

**HAADS/Track 

2.5 

60 

7 

.25 

.1 

- 

- 


LOS 

**HAADS = USAF - High Altitude Attitude Determination System 
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7.2 NOISE MODELS 


The noise model assumed for these analyses identifies the predominant 
noise sources and is based on extrapolating the test data for the 128 x 128 
chip to a 400 X 400 chip design, as presented in Section 6. The resulting 
noise models are tabulated in Table 7-4. 


ACQUISITION MODE 


• Full Frame Scan 

• M y 14 Star Acquisition 

« 2 are-s/s Initial Angular Velocity 

• 0.5 arc-s Error After Acquisition 

• 60 s Acquisition Time 

• No Angular Rate Measurement Available 


TRACK MODE 


• M 7 to M y 14 Stars 

« 0.01 arc-s Accuracy 

• .0028 arc-s Stability 

• 20 Hz Frame Rate 

• 0.2 arc-s/s Maximum Line-of-Sight Rate 

• 5 Star Tracks Simultaneously 


ASTROMETRY 

• M 10 to 17 (and Dimmer ~ Goal of M = 21) 

• Integration Time 50 MS to 600 s v 

• Separate Stars as Close as 0.1 to 1 arc-s 

• Magnitude Accuracy 1% for M = 17 

• Dump All Data in Minimum Time - Typically 2.6 s 

Figure 7-1. Baseline Requirements 


Table 7-4 
NOISE SUMMARY 


Source 

Noise, NEC* 

1. Thermal Noise 

N t = 2.0 t r ' 1/2 . 

t = total read time 
for 1 pixel 
(64 samples) 

2. Dark Current Shot Noise 

N rt , = 51 t, 1/2 
ds i 

t i = time from inject 
to sample 

3. Dark Current Pattern Noise 

N . = 350 t. 
dp i 

t. = time from inject 
to sample 

4. Leakage Current Shot Noise 

"is ■ 25 *r 1/2 

t = time from reset 
to sample 

5. Fixed Pattern Noise 

N f = 2.2xl0 4 

Can be suppressed by 
double-read technique 
to negligible levels 

6. KTC Noise 

N ktc = 5 

Suppressed by cor- 
related dpuble sampl- 
ing and multiple read 
techniques from 
1800 NEC 


*NEC - Noise equivalent carriers, rms, after on-chip differencing, before 
frame-to-frame differencing. 
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7.3 ACQUISITION MODE 

The purpose of the Acquisition Mode is to. reduce the initial state of 
image position uncertainty and velocity to a state which is compatible with 
handover to the Tracking Mode. The Acquisition Mode may also be used to ac- 
quire a second guide star after a slewing command has caused loss of track 
of a first guide star. In this case, a reasonably rapid acquisition will 
be required. The handover state will be achieved by reduction of the errors 
through interaction of the guide sensor, the inertial reference, and the atti- 
tude control system. 

The initial state will depend on characteristics of the ACS and inertial 
reference which are not defined. However, it is assumed that significant 
errors in both velocity and position will be present. 

worst case is assumed on the status of the attitude control at the 
beginning of acquisition, whereby the ACS is started up after a period of 
dormancy and the inertial attitude reference is degraded to the point where 
a full frame acquisition scan is needed to assure acquisition of a suitable 
target star. If the initial attitude is random, there may be no bright star 
in the viewfield. 

During the period of dormancy, the vehicle may have accumulated a signif- 
icant angular velocity. It may be desirable to perform an acquisition before 
stabilizing the vehicle with the aid of an inertial rate reference. These 
assumptions lead to the fol.lowing set of initial acquisition requirements: 

• Full frame scan 

e M y 14 star acquisition 

0 Initial angular velocity - 2 arc-s/s 

0 Final acquisition' error 0.5 arc-s after 60 s 

0 No rate measurement available 

These requirements constitute a worst case situation, but will be used 
to demonstrate how the system can be reactivated after a period of dormancy 
or possibly after partial failure of the inertial reference. 

7.3.1 Effect of Image Velocity 

The integrated signal charge per pixel is inversely proportional to the 
image velocity: 

_ 9 N 

N $p - _p_a ^ where (i) 

2 0 

Qp = pixel subtense angle (0.5 arc-s) 

• 5 

N $ = signal charge generation rate (1.4 x 10 e/s, = 14), 

0 = 1 ine-of-sight rate (arc-s/s) 
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and the factor of 2 allows for the image being split on 2 pixels. N™ may 
vary with the orientation of the image track relative to the pixel structure. 
The maximum per pixel signal charge will be achieved when the image has moved 
a distance equal to two pixel subtenses, or in an integration time 


t 


e 



( 2 ) 


7.3.2 Effect of Position Error 

The initial scan area must be large enough to assure acquisition of 
a detectable star. The density of stars brighter than M v 14 is at least 100 per 
square degree (galactic pole). Each 400 x 400 pixel sensor subtends (400 x 
0.5)^ = 40,000 square arc-s, and six sensors will see an average of 

40,000 x 6 x 100 ^ 2 stars @ M 14 
(3600 T 

To assure acquisition in random orientation, each sensor should therefore 
execute a full frame scan. 

7.3.3 Frame Time 

Consider a full frame scan including a double read operation to cancel 
fixed pattern noise. If we were to read out each 4x4 group as one pixel, 
then the full frame time would be 

t f = 2 x 100 x 100 t r + t g , where 

t = per pixel read time 

t = added time to expose before start of readout (3) 

and two read operations are executed per pixel. For the more critical case 
of high image velocity, the exposure time would be relatively short so that 
t could be neglected in (3). 

The noise model of Table 7-4 provides noise estimates in terms of per 
pixel read time, t , and time since injection, t^. These estimates are based 
on single site readout and must be adjusted to accommodate the new condition 
of 4 x 4 group readout. The corrections are 

• thermal noise x 2 increase 

• dark current shot noise x 4 increase 

• dark current pattern noise x 1 increase 

• leakage shot noise x 2 increase 

and result from the fact that a summation of four rows, each consisting of 
four pixels, would be taken. Thus the signal suiruwould contain (4) ' =2 times 

as much thermal and leakage shot noise and (16) ' =4 times as much dark 

current shot noise as an individual pixel. It is more difficult to assess 
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the multiplier for the dark current pattern noise since it is not, in general, 
randomly distributed pixel-to-pixel . High spatial frequency pattern noise 
would be averaged while low spatial frequency components would be less af- 
fected. A multiplier of unity was assumed for this analysis. The reset noise 
and fixed pattern noise components are assumed to be made negligible and are 
disregarded in this analysis. 

In a frame comprising two read operations per pixel, two thermal and 
two leakage noise samples sum incoherently, but only one sample of dark current 
pattern noise and shot noise is obtained, since the dark charge is injected 
before the second read operation; The total per pixel noise (16 sites) is: 

N t 2 = + (25) V U2 J + ( 51 x 4 h m J + ( 36 ° ‘i) 2 

Substituting tf from (3) and neglecting t , 

N t 2 = 64 x 10 4 /t f + 5000 t f / (2 x 10 4 )+ 4.16 x 10 4 t i + 1.23 x 10 5 t. Z (4) • 

Equation (4) expresses per pixel noise in terms of frame time and time 
since injection. If a line is read out, then injected and re-read, the time 
since injection will equal the frame time. 

From (1) 


sp 


= 46 N s /20 = 1.4 x 10 s t f /9 and 


(S/N)' 


■ 2 x 10 


10 t f 2 /e 2 


64 x 10 4 /t f + 4 . 16 x 10 4 t f + 1.23 x 10 5 t f 2 


S/N 

6 


10 6 t f 3 

32 + 2.08 t f 2 + 6.15 t f 3 


(5) 

( 6 ) 

(7) 


This expression has no maximum but increases with tf. From (2),' however, 
the maximum per pixel signal is reached when the image has traversed two pixel 
subtenses. In this case 


tf = 4/e, for maximum signal (8) 

For 0=2 arc-s/s, tf = 2 s, in which case S/N = 149 for a M 14 star. 

The read rate for this condition would be 2 s/2 x 10 4 = 100 ys. Alternatively, 
one could operate at a fixed frame time of, say, one second. In this case 
the image would generate a track of length 

0 

7 ]— — = 0/2 subgroups 

6 P 
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Thus if e were high, the per pixel S/N would be low according to (7), but 
the number of opportunities for detection would be large. Under these con- 
ditions, it is likely that reliable detection could be achieved at angular 
velocities up to 8 arc-s/s, where S/N = 20. There would be 4 exposed groups 
in a frame. 

The first scan would locate the target within a few groups, depending 
on the angular velocity. The next scan could be restricted to a smaller area, 
resulting in a shorter frame time and a more accurate position measurement. 
Depending on the characteristics of the ACS, hand over to track mode would 
be achieved in a few seconds. 

7.4 TRACK MODE 


For the track mode it is assumed that the system must be capable of han- 
dling a visual magnitude range of M v 7 to M v 14 and provide 0.01 arc-s accuracy 
and 0.0028 arc-s stability at a frame rate of at least 20 Hz. The maximum 
line-of-sight rate is unspecified but, for the sake of discussion, a LOS rate 
of 0.2 arc-s/s will be assumed, consistent with the rate specified by Stanton, 
et al.( 7 ) for solar system tracking. Capability must exist for tracking at 
least five stars either simultaneously or sequentially. Differences in the 
time of star centroid location must be accounted for if sequential tracking 
is used. The five star- images will normally be expected to fall on five dif- 
ferent arrays in order to avoid the penalty of increased readout bandwidth; 
however, should several star images fall on one array, the capability would 
exist for simultaneous tracking qf five stars with concomitant bandwidth in- 
crease. 


7.4.1 Accuracy/Stability 

Tracking accuracy depends upon star visual magnitude, pixel subtense 
angle, pixel-to-pixel variations in sensitivity and linearity, system thermal 
noise, and star signal shot noise. Tracking stability is affected by all of 
the above except sensitivity and linearity variations. It is assumed that 
background variations (offsets) from pixel-to-pixel are reduced to an insignif- 
icant level by the double read method, the noise model is as given in Section 
7.2, and the pixel subtense angle is 0.5 arc-s. 


These effects are analyzed in the Appendix to find the accuracy and 
stability levels attainable. The limiting cases of thermal' noise only, signal 
shot noise only, and linearity variations only, are analyzed for the star 
image centered on a 3 x 3 array and for the image located at the corner of 
the center pixel in a 3 x 3 array. The radius of stability (r § ) and radius 
of accuracy (r^) derived from the calculations are restated below and plotted 
in Figure 7-2. 


Case Al-Star centered, thermal noise only: 

r g = (2.02 x 10 ®)(2.5) v arc-s 

Case A2-Star in corner, thermal noise only: 

o M 

r $ = (2.65 x 10 ') (2.5) v arc-s 


( A30) 


(A33) 
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Figure 7-2. Radius of Stability Vs Star Magnitude 
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Case Bl-Star centered, shot noise only: 

c M v /2 

r s = (1.14 x 10“ ) (2.5) arc-s 

Case B2-Star in corner, shot noise only: 

c M v /2 

r s = (1.23 x 10 ~ d ){ 2.5) arc-s 
Case Cl-Star centered, 1 inear ity variations only: 
r ft = .000828 arc-s 

Case C2-Star in corner,linearity variations only: 


(A37) • 


(A41) 

(A46) 


r A = .000837 arc-s (A49) 

The curves of Figure 7-2 show that the specification of .0028 arc-s 
stability is actually limited to a star of visual magnitude M v < 12 by the 
shot noise on the signal and not by system thermal noise. The linearity varia- 
tions are an insignificant contribution to the radius of accuracy, and so 
the combined thermal noise and signal shot noise allow the .01 arc-s specifica- 
tion to be met by stars of magnitude M y <.14. 

7.4.2 Signal Integration Times 


The required signal integration time depends upon the visual magnitude 
of the star to be tracked, the noise, and the desired S/N ratio for detection 
and tracking; or the noise level required for accurate magnitude determination. 
These parameters are shown in Figure 7-3 in order to highlight the perspective 
values involved for a given application. In Figure 7-3, the times required 
to attain given S/N ratios are plotted against stellar visual magnitude. 

Curve A shows the time required for saturation of the center pixel of a 3 x 3 
pixel group, assuming a centered star image. Curves B, C, and D indicate 
the integration times required to attain S/N ratios of 1000, 100, and 10 respec- 
tively, at temperatures of 0 °C and -25 °C. These curves were generated using 
an approximate formula that includes the effect of system thermal noise, signal 
shot noise, and dark current shot noise. The thermal noise was assumed tg 
be 25 carriers rms and the dark current was assumed to be 2 nA/crrr at 25 °C. 

The S/N values used for the calculations represent the total star signal di- 
vided by the noise on a single pixel, as defined for the centroid calculations. 

It is seen from the curves in Figure 7-3, that system thermal noise 
is of significance only for integration times below about 0.1 second and low 
S/N ratios; otherwise shot noise predominates. Reducing array temperature 
below about 0 °C is seen to be useful only at the very long integration times. 

With the required frame rate of 20 Hz, a S/N ratio of 10, and thermal 
noise level of 25 carriers, stars as dim as M v 18 could be acquired and tracked 
with a radius of stability of about 0.2 arc-s. While it is valid to consider 
the reduction of system thermal noise below 25 carriers by several available 
techniques, such as transform readout, background signal storage, etc., the 
benefits would seem to be limited to the short integration times and low S/N 
ratio conditions. 
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The reader is reminded that the results indicated in Figure 7-3 are based 
on approximate numbers and are presented as a rough design tool and not for 
specification. When a detailed error analysis has been accomplished, this chart 
could serve as a guide for selecting operating modes. The possibility exists 
for automatic programming of integration times for optimum system performance. 

Figure 7-3 also shows that integration times larger than 0.05 second 
would result in signal saturation of the center pixel' for ^v7 stars whereas 
M v 17 stars would permit integration times on the order of 500 s before satura- 
tion occurred. For the longer times, array temperature would be a factor. 

7.4.3 Implementation Approaches 

The preceding section discussed the impact of the various performance 
and noise considerations on the system computations required. This section 
provides further detail on computation methods which may be implemented for 
the track mode. The selected method will depend on the particular detail re- 
quirements, but for an application where minimum noise is essential, the recom- 
mended method would be to use a 4 x 4 submatrix array with four parallel ampli- 
fiers as an output to reduce data rates and provide maximum flexibility for 
tracking multiple stars. The basic scheme would use correlated double sampling 
at an optimum sample rate selected to minimize leakage, KTC, and Johnson noise. 
Sixty-four samples would be accumulated for a given read time using the NDRO 
capability of the CIO to enhance the signal-to-noise ratio. Digital double 
correlation would be used to reduce pattern noise by 60 db. A chip temperature 
of 0 °C appears adequate for tracking up through M v 14 stars, but for astrometry, 
further cooling is indicated. To determine if the stability goal of .0028 arc-s 
for M v 14 stars could be met, further analysis is needed. It appears from the 
preceding analyses that the system would require a significant reduction of 
the 20 Hz frame data rate in order to permit sufficient integration time to 
achieve the necessary S/N ratio (*1000) to meet the stated stability goal. 

The proposed mechanization can handle five stars, but with some further band- 
width penalty over a single star, unless further parallel outputting is used. 

A more detailed analysis and design effort is required to establish the specific 
mechanization. 

One of the key features of the CID is the ability to randomly access 
individual pixels. The random access can be achieved by scanning registers 
on the chip or by using decoders to address a specific pixel in a shorter time 
period. The choice of one approach over the other again depends on the specific 
implementation requirements. 

7.5 ASTROMETRY MODE 


The astrometry mode is assumed to require a readout of stars of between 
M v 10 and M v 17 with a desirability of reading out up to M v 21. Integration time 
requirements vary from a minimum of 50 ms to as long as 600 s, the maximum 
integration time required for a magnitude 21 star. 

The system should be capable of separating stars spaced as close as 
0.1 arc-s, and of determining magnitude to an accuracy of 1 percent for a 
M y 17 star. Finally it should be capable of reading out all the pixels in 
2.6 s, based on FGS requirements. 
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7.5.1 Variable Magnitude 

The CID appears to be able to handle all of the star magnitudes indicated 
in the requirement. The dimmer stars may require additional cooling in order 
to avoid dark current saturation and to achieve accuracies for the magnitude 
required. The CID is unique in the fact that it is possible to allow the sys- 
tem to integrate while reading. Even if individual pixels saturate, those 
which are not saturated can be read. 

7.5.2 Separation of Stars 

The requirement indicates a desirability of separating stars as close 
as 0.1 arc-s. Based on the current design with a pixel of a half arc-s and 
typical blur circle approximately the same size, the presence of two stars 
one-tenth of an arc-s apart would cause a change in the size of the blur cir- 
cle. It appears to be possible to measure this difference in the order of 
20 percent' and hence detect star separation in the order of 0.1 arc-s as is cur- 
rently done with. photographic emulsions. 

7.5.3 Magnitude Accuracy * 

The requirement for 1% magnitude accuracy requires that the signal to 
noise be at least 100 after a multiple reading and that the sensor responsivity 
variations be 1% or less. This performance appears to be within the capabilities 
of the system. Figure, 7-3 indicates that adequate signal to noise ratios appear 
to exist for star magnitude up to M 18. Pixel-to-pixel variations of <1 percent 
have been measured on representative arrays. If necessary, calibration could 
be used for 'the pixels of interest. All pixels could be stored in memory but 
this would take more memory than desired. In this case, selected subsets of 
the array could be useful if the vehicle maneuvered such that the stars of 
interest appear in these calibrated pixels. 

7.5.4 Complete Array Dump 

•The requirement for reading all of. the pixels in as short a time as pos- 
sible has several ramifications. In the JPL report, ( 8 7 2.6 s was used at a time 
to dump the charge in all pixels. If this were used for the 400 x 400 CID 
(.06 ms/pixel), it would increase the read rate over that for the track mode 
(.003 ms),. The temporal noise would then be as large as 100 carriers for this 
rapid readout and hence loss of signal to noise would occur. An alternative 
is to read at the slower rate. It would take as long as 50 s to read all the 
pixels. If selected areas were read out, the time would be a function of the 
number of pixels involved. Several options exist as to how to read all of the 
data in the array. 

7.5.5 Implementation 

Based on the above discussion it appears that in the astrometry mode 
many of the techniques that have been discussed for the track mode would also 
be used. However, the main difference is that the readout rate requirement 
is open and increased integration times are acceptable. It is possible with 
the CID to continue integration after some pixels have saturated while si- 
multaneously reading out other pixels of interest. This feature provides 
a major asset for the astrometry mode because of the ability to improve 
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dynamic range. Specifically, the pixels may be read out at a rate correspond- 
ing to a minimum noise rate. Repetitive reading will enhance accuracy and 
permit linearity measurements. The ability to read pixels at random allows 
one to keep track of background current, dark current, and variations in the 
star signals over the total integration time. This provides a unique CID 
capability. 

The exact method of reading out the data to telemetry is a matter of 
the particular interface details. The JPL report (9) shows the use of two 
line buffers in order to store data prior to transmitting over the telemetry 
links. A similar approach would be anticipated for the CID but because of 
the ability to separate readout from integration, the readout times are not 
critical . 

7.6 COMPUTER IMPLICATIONS 

A brief review of the computational implications was made in order to 
assess the major areas of concern in implementing the computational algorithms 
In general the primary forcing functions for the star trackers are in the 
readout rate areas and readout accuracies. The general problem of computer 
memory to store the required algorithms and to store data appears to be a 
relatively minor problem in light of the large memory capacity of existing 
technology. The computations required are relatively simple and have little 
impact on the throughput capacity of current microcomputers. 

Typical numbers depend so much on the application that it is difficult 
to generalize. However, the JPL study of the Fine Guidance Sensor indicated 
on the order of 1,000 words of RAM and 1,000 words of PROM would be adequate. 

A similar study for a different application by Rockwell on a High Altitude 
Attitude Control Demonstration System (HAADS) utilized a fairly complex fit 
pattern threshold calculation and stored a great deal of data in order to 
minimize this dark current bias correction factor. In this case a random 
access memory of 4,000 words and a PROM memory of 4,000 words was required. 

A brief look of the Fine Guidance System and the computations that we would 
anticipate for most applications, indicate that a maximum of 4,000 words of 
RAM and 4,000 words of PROM memory are all that would normally be required 
except where excessive amounts of calibration of individual pixel data are 
desired. For example, if in the case of a 500 x 500 pixel array, all pixels 
were to be calibrated with ah 8 bit word this would require 250,000 words, 

2 x 10° bits. Using today's technology, 64,000 bits on a single chip are 
reasonable ROM memory. In this case, the amount of memory would take about 
30 chips which is a relatively small hardware penalty. 

The accuracy requirements tend to fall into two categories. For most 
computations, accuracies in the order of 8 bits are ample. However, for ab- 
solute value calculations, 16 to 18 bits are often required. For example, 
if a quantization of .001 arc-s is needed to minimize errors in stability 
calculations, then for a pixel subtending one-half arc-s in an array of 500 
pixels, 18 bits is required. It thus appears that 18-bit accuracy would be 
needed in many cases. This is consistent with the JPL approach which pro- 
vides for 18 bit input and output data, needed for the precise centroid cal- 
culations. 
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In summary, the total computational load impact is relatively small 
and can be handled with minimal impact on the hardware. It is necessary to 
keep track of the A/D rates and the A/D accuracy requirements in order to 
minimize system cost and complexity. 
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Section 8 

8.0 CONCLUSIONS AND RECOMMENDATIONS 


The test results, coupled with recent advances in silicon device proces- 
sing capabilities, clearly indicate that the CID has great potential as a stel- 
lar tracking device in applications with performance requirements such as those 
for the Space Telescope Fine Guidance Sensor. 

Projected performance parameters of a large area CID array designed for 
star tracking are: 

• 20 pm x 20 ym square sensing site; 400 x 400 array 

9 Noise <25 carriers/pixel 

• Quantum Yield >40%, from 0.4 to 1.0 ym 

• Modest cooling requirements, on the order of 0°C 

o Low potential cost due to structural simplicity 

• Simple peripheral circuitry 

Random access and on-chip processing capabilities are potentially useful for 
further system simplification. 

Further effort is indicated to verify the projected CID performance. 

This includes: 

1. Design, fabrication and test of a full 400 x 400 array, including 
the necessary system definition, low-noise readout electronics, 
high accuracy data acquisition circuitry, software, etc., to fully 
evaluate performance in the star tracker mode. 

2. Fabrication and test of a breadboard star tracker module using 
the array developed in 1, above, or as an interim .approach, using 
an existing array design. Using appropriate star image simulation, 
the breadboard module should be evaluated against previously de- 
fined system functional requirements. 
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Section 1 0 


10.0 APPENDIX 

10.1 Calculation of Tracking Accuracy and Stability 

This appendix presents the calculation of the radius of stability and 
accuracy for the limiting cases according to the following specifications: 

1. A 3 x 3 pixel matrix is used for centroid calculation. 

2. The star image is gaussian with a spread of ex = one pixel width. 

S p 

3. The effective aperture area of the lens is A = 32,000 cm . 

4. The amount of stellar radiation reaching the lens for a M 0 

star, I Q = 3.57 x 10 6 photons/cm2-s. v 

5. The integration time is t '= .05 s. 

6. The combined lens and imager quantum efficiency over- the portion 
of the spectrum of interest is Q‘ = 0.5 carriers/photon. 

7. One pixel width subtends an angle of .5 arc-s. 

8. The level of random temporal noise is 25 carriers. 

9. The linearity of response variations are 0.5%. 

The centroid calculation is performed on the 3 x 3 pixel matrix shown 
in Figure A-l. 



-1 0 1 

/ 


Figure A-l. 3x3 Pixel Sub-array 
The equations used to calculate the centroid position are: 


+1 
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where V.. is the signal at pixel location (i,j) and 
^ 3 

+1 +1 

T = £ E V,, (A3) 

i-1 j-1 

is the total signal falling on the 3x3 array. 

For large signal-to-noise ratios, the variation in the X direction 
is 

dX = E E (3X/3V . .)dV . . (A4) 

i j • 

where (3X/3V-.) * (3/aV..) (1/T)(E E i V..) 

ij (J [ -j j '0 ^ 

= ( 3/3 v. .) (E I j V..)/(E E V..) 

J L 1 j J 1 j 

= (j E E V - E E j V )/& 
i j J i J V 

= (jT - E E j V )/T 2 
i J J 

= &L (A5) 

The best case for accuracy and stability is when the star is centered 
in the matrix and the spread of light is as shown in Figure A-2. The worst 
case, assuming the 3x3 block can be moved in single-pixel increments, is 
when the star is in the corner of the center pixel and the spread of light 
is as shown in Figure A-3. 

1 

<1 0 

-1 

-i je i 

/ 

Figure A-2. Light Distribution - Star in Center 
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10.1.1 Centered Star 

For the case of the centered star, the variance in the X direction is 

dX 2 = (1/T 2 ) E E (j-X) 2 o N 2 (i,j) (A6) 

* j 

2 

where (i,j) is the noise on the (i,j) pixel. In this case 

X = 0 .and Y = 0 (A7) 

The variance in X is 

dX 2 = (1/T 2 ) / 2 [( -1 ) 2 o n 2 (C) + (0) 2 o N 2 (B) + (1)? a„ 2 (C)] 

+ [(-l) 2 a/(B) + (0) 2 a/ (A) + (l) 2 a N 2 (B)]| ■ 

= (1/T 2 ) 1 4 o n 2 (.C) + 2 o N 2 (B)| (AS) 

The star's intensity profile is the standard two-dimensjonal normal 
distribution, 

V(X, Y) = (V Q /27ra s ) exp [(X 2 + Y 2 )/2 a g 2 ] (A9) 

where V Q is the total star radiation as measured by the imager, 

V Q = (I Q A Q t)/(2 .5) v (A10) 

When the proper values are inserted, this becomes 

- M 

V Q = (3.57 x 10 b ) (3.20 x 10 4 )( .5) ( .05)/(2.5) v 

q M 

= 2.86 x 10 /(2. 5) v (All) 
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where M v , the visual magnitude, has been left as the independant variable. 
The star signal which falls on the 3x3 array is 


T = f 1 ’ 5 f 1,5 V(X,Y) dXdY 

J- 1.5 J- 1.5 


= fl.5 r 1.5 (y /2 7rcf j exp (xV)/2a^ dXdY 
j- 1.5 •'-1 . 5 0 s L M 

= (V 0 /27ra $ ) r 1 - 5 exp X 2 /2 o\ fl.5 exp y 2 /2 dYdX 
J -1.5 L J J_ 15 L SJ 

= V /(2ira ) 1/2 1 f 1,5 exp|"x 2 /2 a 2 j f*l/( 2tt) 1/2 1 f 1 ' 5 expfY 2 /2a 2 (dY/a )l|d 

L° s J j -i.5 L s J[L _p-i . 5 L J 

or with the substitution Z = Y/a $ , this becomes 

T = V o (27ra s ) 1/2 f 1-5 exp X 2 /2a s 2 | l/(2m) 1/2 exp Z 2 /2 dzjdX (A1 

L ^ -1.5 L JL -1.5 . L J 


= V o /{2 7ra s ) 1/2 J 1 * 5 exp X 2 /2a g 2 erf^j - erf dX 

The same operation can be used for integrating with respect to X. The result 
is 


T = V o ^erf(1.5/a s ) - erf -(-1.5/a s )J‘ (Al< 

For a s = jl pixel width this can be evaluated with the help of Table A-l. The 
result is 


■v 


433 + .433 


M ■' = .750 V c 

The variaiance in the x direction is 


dX 2 = 4a N 2 (C) + 2a N 2 (B)]/(.750 V Q ] 
*“ J 


The radius of stability is 


r s = (dX 2 + dY 2 ) 1/2 


(A17) 
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Table A-l 

TABLE OF THE ERROR FUNCTION 


er'fx = 




X 

/o 


exp(-y2/2)dy 


X 

erf x 

X- 

erf x 

0.05 

0.01994 

1.55 

0.43943 

0.10 

0.03983 

1.60 

0.44520 

0.15 

0.05962 

1.65 

0.45053 

0.20 

0.07926 

1.70 

0.45543 

0.25 

0.08971 

1.75 

0.45994 

0^30 

0.11791 

1.80 

0.46407 

0.35 

0.13683 

1.85 

0.46784 

0.40 

0.15542 

1.90 

0.47128 

0.45 

0.17364 

1.95 

0.47441 

0.50 

0.19146 

2.00 

0.47725 

0.55 

0.20884 

2.05 

0.47982 

0.60 

0.22575 

2.10 

0.48214 

0.65 

0.24215 

2.15 

0.48422 

0.70 

0.25804 

2.20 

0.48610 

0.75 

0.27337' 

. 2.25 

0.48778 

0.80 

0.28814 

2.30 

0.48928 

0.85 

0.30234 

2.35 

0.49061 

0.90 

0.31594 

2.40 

0.49180 

0.95 

0.32894 

2.45 

0.49286 

1.00 

0.34134 

2.50 

0.49379 

1.05 

0.35314 

2.55 

0.49461 

1.10 

0.36433 

2.60 

0.49534 

1.15 

0.37493 

2.65 

0.49597 

1.20 

0.38493 

2.70 

0.49653 

1.25 

0.39435 

2.75 

0.49702 

1.30 

0.40320 

2.80 

0.49744- 

1.35 

0.41149 

2.85 

0.49781 

1.40 

0.41924 

2.90 

0.49813 

1.45 

0.42647 

2.95 

0.49841 

1.50 

0.43319 

3.00 

0.49865 
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and since dX 2 --dY 2 at the center, 


r $ = (2dX 2 ) 1/2 


Sa N 2(c) + 4a N 2 {B) 


1/2 


/( . 750 V ) 


(A18) 


10.1.2 Star in Corner 


For the star in the corner of the- center pixel, equation (A6) still holds : 
dX 2 = (l/T 2 )Sp(j-x) 2 CT n (i,j) (A19) 


In this case 


X = .5 and Y = .5 

The variance in the X direction is 


dX 2 = {1/T 2 )<2 


(-1.5) 2 o n 2 (E) + (-.5) 2 o n 2 (D) + (.5) 2 o n 2 (D) 


(-1.5) 2 a N 2 (F) + (.-.5) 2 c> n 2 (E) + (,5) 2 a N 2 (E) 


» (1/T 2 )|2.25o n 2 (F) + 5o n 2 (E) + a^D) J 
The star signal which falls on the 3x3 array in -this case is 


} 


T = J 1 J 1 V(X,Y) dXdY 

= V 0 [erf(l/a s ) - erf(-2/a s )] 2 
With the use of Table A-l, this evaluates to 

T = V q [.341 + .477 
= 0.669 V 

o 

The variance in the X direction is 


dX 2 = [2.25 <j n ( F) + 5a N (E) + ° N (F)] /{.669 V Q ) 2 
The radius of stability js 
r s = ( 2dX 2 ) 1/2 

= [ 4>5a N (F) + 10 °N (E) + 2 ct n (D)] 1/2 /(-699 V q ) 


(A20) 

(A21) 

(A22) 

(A23) 

(A24) 

(A25) 

(A26) 
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10.1.3 Case A 

This case calculates the radius -of stability when there is only random, 
uncorrelated, thermal noise, and assumes that there is.no shot noise.' There- 
fore 

ct N (A) = ct N (B) = ’ = a N (F) = a N (A27) 

The radius of stability for a star that is centered. Case Al, is from equa- 
tions A18 and A27 

r S * (12 " n 2 ) 1/2 /(-750V o ) 

= 4.62 (o N /V 0 ) (AZ8) 

With the noise estimated to be 25 carriers and V as given by equation All, 
the radius of stability is 0 

M 

r<~ = (4.62) (25) (2.5) v /2.80 x 10 9 
o M 

= (4.04 x 10 °)(2.5) v pixel widths (A29) 

. In arc-s this is 

o M 

r s = (2.02 x 10"°) 2.5 v arc-s (A30) 

This is plotted in Figure 7-2, Curve -Al. 

The radius' of stability for a star centered at the corner of the center, 
pixel. Case A2, is from equations A26 and A27 

r. = (16.5a 2 ) 1/2 /(. 669V ) 

•-> N o 

= 6.07 (a M /V o ) ( A31) 

When the proper values are inserted this gives 

M o 

r<- = (6.07) (25) (2.5) v /2.86 x 10 9 

o 

= (5.31 x 10 ) (2 . 5) v pixel- widths (A32) 

This is equivalent to 

8 

r s = (2.65 x 10~°)(2.5) v arc-s (A33) 

= 1.32 times Case Al 
This is plotted in Figure 7-2, Curve A2 
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10.1.4 Case B 


The other limiting case assumes that there is no thermal noise, only 
shot noise on the star signal. In order to find the radius of stability 
for this case, each individual noise level must be calculated from 


ycj) 


V(i,j) 


1/2 


(A34) 


When the star is centered. Case Bl, V(B) and V(C) must be calculated. They 
are 


V(B) = 


r.5 rl.5 V ( XjY ) dXdY 
J -.5 J .5 


= V o [erf(.5o $ )-erf(-.5a s )] [erf (1.5/o s )-erf ( . 5/a s )J 



.191 + . 19lj [.433 - . 19l] 


= .0924 V 0 


(A35) 


and 


1.5 1.5 

V(C) = f f 
.5 J .5 


V(X, Y) dXdY 


= V Q [erf(1.5/a s ) - erf(.5/a $ )] 2 
= V o [.433 - . 19 lj 2 
= .0586 V Q 

'The radius of stability is from equations A18, A34, A35, and A36 


r s = [8 (.0586V q ) + 4 ( .0924V q )] 1/2 /( .750 V Q ) 
= 1.22/(V o ) ly/2 

= [(1.22) (2.5) v /2 ]/(2.86x!0 9 ) 1/2 

M /2 

- (2.28 x 10 -5 ) (2.5) v pixel widths 
, _ 5 M v /2 

= (1.14 x 10 °)(2.5) arc-s 

This result is plotted in Figure. 7-2, Curve Bl. 


(A36) 


(A37) 


When the center of the star falls in the corner of the -3 x. 3 array, 
Case B2, V(D), V(E),. and V(F) must be calculated. They are 

V { D ) = f 1 f 1 V(XY)dXdY 

J 0 Jo 
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and 


- V 0 [erf(l/a s )] 2 

= V [. 34l| 2 
o[ 1 

= .116 V Q ' (A38) 


v ( E) = J 1 j* 2 V ( XY ) dXdY 

= V 0 [erf(l/a s )] |erf(2/a s ) - erf(l/<^)] 

= V Q [.34l] [.477 - . 34l] 

= .0464 V 0 (A39) 


CHIU p p 

V(F) = / 1 / 1 V ( XY ) dXdY 

= V o [erf(2/a $ ) -erf{l/a s )] 2 

= V Q [.477 - .341 

= .0185 V (A40) 

o- 

The radius of stability is, from equations A26, A34, A38, A39, and A40 
r s = [4.5(.0185 V Q ) + 10(.0464 V Q ) + 2 (.116 V Q )j 1/2 /( .699 V Q ) 

= 1.32/(V q ) 1/2 

= [(1.32)(2.5) Mv/ ]/(2.86 x 10 9 ) 1/2 
-5 M v /2 

-= (2.47 x 10 )(2.5) v pixel widths 

-5 V 2 

= (1.23 x 10 °) (2.5) v arc-s 


= 1.08 times the result of Case B1 - (A41) 

This result is plotted in Figure 7-2, Curve B2. 

10.1.5. Case C 

For this case the radius of accuracy is calculated. In addition to the 
random noise results of Case A and Case B, there will be a fixed pattern noise 
due to (1) random effects such as threshold variations and (2) linearity of 
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response variations. The double read method is expected to cancel the random 
(but positionally stable) fixed pattern noise to the extent where is insignifi- 
cant. The' linearity of response variations are estimated to be . 5 % or 

a = .005 ( A42) 


and so the noise is 

a n (XY) = a V(XY) 

A 


( A43 ) 


When the star is centered. Case Cl, the radius of accuracy is, from 
equation A18, 


r ft = [ 80 n 2 (C) + 4 a n 2 (B)] 1/2 /(.750 V Q ) (A44) 

With the use of equation (A43) this becomes 

r A = [s a 2 V{Cy 2 + 4a 2 V^Bj 2 ] 1/2 /(.750 V c ) (A45) 


and from equations A35 and A36, this becomes 

r A - = [8a 2 (.0586 V q ) 2 + 4 a 2 (.0924 V q ) 2 ] 1/2 /( .750 V Q ) 

= a |. 33 1] 

= .005 [.33l] 

= .00166 piixel widths 

= .000828 arc-s (A46) 

This result is plotted in Figure 7-2, Curve Cl 


When the star is centered in the corner of the center pixel, the radius 
of accuracy is, from equation A26, 


4.5a n 2 (F) + 10a n 2 (E) + 2a n 2 (D)] /( .669 V Q ) 


(A47) 


With the use of equation (A43) this is 

r A = [ 4 * 5 a2 + 10a2 T ^ 2 + 2a2 VTDj 2 ]/( -66 9 V Q ) (A48) 

When this result is combined with equations A38, A39, and A40, it gives 

r n = [4.5 a 2 ( .0185 V } 2 + 10 a 2 (.0464 V ) 2 + 2 a 2 (.116 V ) 2 l 1/2 /{.669 V ) 

A l 0 0 0 0 
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= a (.335) 

= ( .005) ( .335) 

= .00167 pixel widths 
= .000837 sec 

( A49) 

= 1.01 times the result of Case Cl 
This result is plotted in Figure 7-2, Curve C2. 
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